Non-locality in the nucleon-nucleon interaction and nuclear matter saturation 
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We study the possible relationship between the saturation properties of nuclear matter and the 
inclusion of non-locality in the nucleon-nucleon interaction. To this purpose we compute the satura- 
tion curve of nuclear matter within the Bethe-Brueckner-Goldstone theory using a recently proposed 
realistic non-local potential, and compare it with the corresponding curves obtained with a purely 
local realistic interaction (Argonne vig) and the most recent version of the one-boson exchange 
potential (CD Bonn). We find that the inclusion of non-locality in the two-nucleon bare interaction 
strongly affects saturation, but it is unable to provide a consistent description of few-body nuclear 
systems and nuclear matter. 

PACS numbers: 21.65.+f, 24.10.Cn, 21.30.-x 



I. INTRODUCTION 

The effective nucleon-nucleon (NN) interaction is ex- 
pected, on physical grounds, to be intrinsically non-local. 
The quark structure of nucleons implies that the po- 
tential between two interacting nucleons cannot be ex- 
pressed only in terms of the center of mass degree of free- 
dom, at least for distances smaller than twice the nucleon 
radius. Even a description of the NN interaction in terms 
of one-boson exchange processes introduces some degree 
of non-locality []]. For practical purposes, however, it 
is convenient to have a local (energy independent) NN 
interaction, to be used in numerical applications. The 
modeling of the NN interaction by a local potential, with 
possible guidance from meson exchange processes, has 
been extremely successful in reproducing the NN exper- 
imental phase shifts and deuteron properties 0, ■ This 
approach has developed for many years and has reached 
a high degree of sophistication. 

It has been established by now that these local NN 
potentials, which are essentially phase equivalent, give 
slightly different results for the three- and four-nucleon 
systems. Furthermore, in general they underestimate the 
binding energy and do not give the correct values for 
some polarization observables. This drawback has been 
usually overcome by introducing three-body forces, which 
can be again justified by the non-elementary nature of 
the nucleon. Actually three-body forces can be gener- 
ated also by meson-exchange processes of higher order, 
where the nucleon is excited in the intermediate states of 
the interaction processes or meson-meson couplings are 
introduced p|. Both non- locality and three-body forces 
can be considered to have a common origin, the inner 
structure of the nucleon. However, one cannot fully re- 
duce one to another, since non-locality has definite effects 
already at the two-body level, e.g. a non-local potential 
cannot be equivalent, in general, to an energy indepen- 
dent local potential. Furthermore, it is not clear to what 
extent the off-shell behavior of a non-local potential, i.e. 
the off shell T-matrix, can be simulated by the presence 
of three-body forces. 

Another question where non-locality could play a role 



is the mechanism of saturation in nuclear matter. In 
fact, it turns out that in order to get saturation with 
a local (energy independent) interaction, it is essential 
to include a strong repulsive core, as also suggested by 
the behavior of the NN phase shifts. Again, three-body 
forces are necessary to tune the saturation point close to 
the empirical one. It is surely conceivable that a certain 
degree of non-locality could simulate the strong repulsion 
at short distance, which is one of the main features of the 
NN interaction, and maybe reduce the need of three-body 
forces. 

It appears, therefore, of great interest to investigate to 
what extent non-locality alone could enable to get satura- 
tion, possibly close to the empirical one, and at the same 
time to describe three- and four-body systems correctly. 

Recenth/jin a series of papers by Doleschall and collab- 
orators [jj 0, > a new realistic non-local NN potential 
was proposed, which is able to describe accurately both 
the NN phase shifts and the properties of three-nucleon 
systems ( He and 3 H), i.e. their binding energies and 
radii. Also the alpha particle appears to be reasonably 
well described |8j . The need of three-body forces in this 
case seems to be absent, or completely simulated by the 
non-locality. 

This NN interaction introduces non-locality at short 
distance (r < 3 fm _1 ) for the s-wave channels, and possi- 
bly also for the p- waves ones Pj . The rest of the potential 
is taken as the Argonne vig |2| , while the non-local part 
is purely phenomenological. 

In this paper we address the issue of the saturation 
properties of non-local potentials. To this purpose we 
present nuclear matter calculations within the Bethe- 
Brueckner-Goldstone method extended up to three hole- 
line contributions, which is known to converge well and 
to give an accurate nuclear matter saturation curve. 

As the main representative of non-local NN interac- 
tion we take the potential of ref. Q. We compare the 
corresponding nuclear matter Equation of state (EOS) 
with those obtained starting from the charge-dependent 
(CD) Bonn |9| and the vig potentials. The former can 
be considered the most recent and accurate interaction 
based on the one-boson exchange model and, therefore, 
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it contains some degree of non-locality. The latter is a 
convenient reference local potential for the comparison 
with Doleschall's, which, as we stated above, coincides 
with the vis in the L > 1 channels. Additionally, other 
modern local potentials, like Argonne V14 Q or Paris po- 
tential 0| , give similar results, close enough to vis to be 
considered equivalent for the present purpose 

The study of phase-shift equivalent NN interactions 
and their predictions of nuclear matter saturation prop- 
erties has a long history. For example in refs. [13, [13] 
unitarily equivalent interactions which differ in the short 
range part of the 1 Sq channel were shown to provide 
rather different values of the saturation energy and den- 
sity. More recently, the predictions of modern NN po- 
tentials were studied in refs. [Til ITEl ITtSj. Here the role 
of non-locality in the binding energy of nuclear matter 
was also partially discussed, using the CD Bonn and the 
Nijmegen-I [T^ | potentials, which were shown to provide 
a stronger binding than local potentials such as the vis- 
However, both these potentials underestimate the bind- 
ing energy of three- and four-body systems |19|. at vari- 
ance with Doleschall's interaction [||. Therefore when 
using the latter we can, in principle, expect different re- 
sults also for nuclear matter calculations. 

It is well known that the saturation points obtained 
within the Brueckner approximation for the gap choice 
and different NN interactions tend to lie along a band 
[Tl| , the celebrated "Coester band" , in the energy-density 
plane. The position inside the band for a given NN inter- 
action turns out to be correlated with the d-state proba- 
bility in the deuteron. The high-energy, high-density part 
of the band is usually associated with low values of this 
probability. However, as shown in ref. |18| . this trend 
does not hold any more if three hole-line contributions 
are included, as we do in the present paper. In this case 
the results are quite insensitive to the d-wave probability 
value and actually close to each other, at least for lo- 
cal potentials or quasi-local potentials like the Paris one. 
The same conclusion holds true if the continuous choice 
for the single particle potential is used in the Brueckner 
approximation 

The paper is organized in the following way. In Sec- 
tion [n] we shortly study the free space properties of the 
Doleschall, CD Bonn and vig potentials. In Section HTT1 
after introducing the BBG formalism, we study the EOS 
of nuclear matter obtained with Doleschall's potential, 
addressing, in particular, the issue of the convergence 
of the hole-line expansion in this case. We then com- 
pare this EOS with those obtained starting from the two 
other potentials. Finally in Section Hvl we draw our con- 
clusions. 



II. LOCAL VS. NON-LOCAL NN POTENTIALS. 

Before studying infinite nuclear matter and in order 
to better understand the role played by the different NN 
interactions in the nuclear matter calculations of next 
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FIG. 1: Diagonal matrix elements of the vis (solid line), 
CD Bonn (dotted) potentials and of DoleshalPs IS potential 
of ref. |g (dot-dashed, labeled as NL), for the np interaction 
in the 1 So (upper panel), 3 Si (middle) and in the coupling 
3 Si- 3 Di (lower panel) channels. The on-shell R(qo,qo', E) R- 
matrix elements are also shown (symbols), for laboratory en- 
ergies ranging from 50 to 350 MeV, in steps of 50 MeV. 



section, it is useful to briefly compare the properties and 
characteristics of the potentials we are going to employ, 
on the basis also of some established results obtained in 
the literature HHEl 

The general matrix elements in the momentum repre- 
sentation of a NN interaction V c (r,r') in a given two- 
body channel c can be written 



^drr^dr'jiikr^, (r, r')j[{k'r') (1) 



where Z, V are the initial and final orbital angular mo- 
menta in the channel, and the corresponding spher- 
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FIG. 2: Left panels: off-diagonal 1 So np matrix elements of the vig, CD Bonn and NL potentials, for three fixed momenta 
go corresponding to laboratory energies Ei a t = IE = 25 (upper panel), 150 (middle) and 350 (lower) MeV. Right panels: 
corresponding half off-shell i?-matrix elements. The solid dot marks the on-shell values of the i?-matrix. 



ical Bessel functions. For a local potential V c (r, r') = 
v c S(r — r')/rr'. In the momentum representation the 
matrix elements Vft,{k, k') of eq. are in any case a 
function of two variables, the initial and final relative 
momenta k and k' . The non-locality is incorporated in a 
non trivial way in the analytical properties of the matrix 
elements of the potentials, as well as of the correspond- 
ing T-Matrix, or of the real i?-Matrix, which is used here 
for later convenience. The latter describes the free space 
NN scattering and is obtained by solving the Lippmann- 
Schwinger equation, which for two nucleons interacting 
with energy E in the center of mass frame can be written 
as: 

Rf l ,(k,k';E) = V l Uk,k , )+ (2) 
r°° m 

+X> J o 9 2 dqV l i"(k,q) ME _ q2 m > Aq,k';E), 

where the symbol V denotes the principal value integral 
and M is the free nucleon mass. 

For realistic potentials the on-shell i?-matrix, 



R(qo, <7o; E) with qo — V ME, is fixed, to a large extent, 
by the fitting of the experimental NN phase shifts. 

The potentials we consider in the prese nt p aper were fit 
to reproduce the Nijmegen phase shifts [2(| for energies 
up to 350 MeV in the laboratory frame, and therefore 
give rise to practically identical on-shell values of the R- 
matrix in this energy range. However, due to their totally 
different analytical structure, a quite different behavior of 
the bare NN matrix elements of eq. QJ can be expected, 
especially when non-local potentials are compared with 
local ones. 

This is illustrated in fig-Q where the bare NN diagonal 
matrix elements and the corresponding on-shell i?-matrix 
for the local vj.8 and the non-local CD Bonn are com- 
pared with those obtained from the non-local potential 
indicated as IS in ref. and hereafter referred to as NL. 
Here the dominant S-wave channels are considered, for 
which the main non-locality has been introduced in ref. 
@- One can see that, despite the expected agreement of 
the various interactions on the i?-matrix, the behavior of 
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the bare NN matrix elements is quite different. At small 
momenta we see that in the 1 So channel the two non-local 
potentials are rather close to each other, while they dif- 
fer significantly from the Vis This is still true in the 3 Si 
channel, where, however, the difference between NL and 
CD Bonn is increased. It can be also seen that at increas- 
ing non-locality, i.e. going from the CD Bonn to the NL 
potential, the deviation form the local potential vis is in- 
creasing. This is in agreement with refs. [Til HE}, where 
it was shown that including non-locality makes the at- 
traction provided by the integral term of eq. weaker, 
thus requiring a less repulsive diagonal matrix element 
of V, in order to obtain the same on-shell i?-matrix. In 
particular it is interesting to note that the diagonal NL 
potential is always rather close to the corresponding on- 
shell i?-matrix elements, especially for the 3 Si channel. 
In the coupling channel 3 Si- 3 .Di we see that the three po- 
tentials are very close to each other at small momenta, 
with NL deviating from the others at intermediate q. 

Considering the large momentum behavior, we observe 
that in all cases the CD Bonn potentials exhibits a much 
larger tail. 



Since the diagonal matrix elements of the bare poten- 
tials are already so different from each other, is not sur- 
prising then that the off-shell behavior of both the bare 
V and the i?-matrix is also different, as illustrated in 
figs. [3 and |2| where the half off-shell matrix elements of 
the bare potentials and of the i?-matrix are shown for 
three choices of the incident energy. The behavior of the 
bare potentials (left panels in the figures) shows essen- 
tially the same features we have already observed for the 
diagonal matrix elements. The curves representing the 
i?-matrix (right panels) are closer to each other, but still 
sensible differences are observed, especially when the re- 
sults corresponding to the vis are compared with the two 
other non-local potentials. 

Finally in fig. 0] we also show the half off-shell V and 
R 3 S-Di matrix elements. Here only the tensor part of 
the NN interaction can contribute, which was shown in 
previous works to be strongly affected by the inclusion 
of non-locality 0, 0, , and which is known to play 
an important role in nuclear structure calculations. We 
see that indeed, the inclusion of non-locality reduces the 
strength of the interaction and, for the NL case, also the 
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range of momentum values where the matrix elements 
are significantly different from zero. 



III. EOS OF NUCLEAR MATTER 

A. The BBG expansion 

The microscopic calculations of nuclear matter Equa- 
tion of State, i.e. the energy per particle as a function of 
density, on the basis of realistic interactions has a long 
history. Theoretical and numerical methods, accurate 
enough for our considerations, are nowadays available. 
We will follow the Bethe-Brueckner-Goldstone scheme, 
which has been proved [Til ] to be reliable in a wide range 
of density. In this expansion the original bare NN in- 
teraction v is systematically replaced by the G-matrix, 
which satisfies the integral equation 



G[p; u>] = v + v 



\kgh) Q (kgh\ 

lj - e(k a ) - e(kb) 



G\p-M 



(3) 



where e(k) is the single particle energy and Q the so- 
called Pauli operator which projects both intermediate 
single particle states above the Fermi surface. The di- 
agrams in terms of the G-matrix are then ordered ac- 
cording to the numbers of hole-lines they contain. The 
introduction of the self-consistent single particle poten- 
tial U (k) is an essential ingredient of the method 



U{k)= (kk'\G{e(k) + e{k'))\kk') 



(4) 



k' <kf 



Since e(k) = h 2 k 2 /2M + U(k), eqs. (EH) include a self- 
consistent procedure which determines U(k). In particu- 
lar we use the continuous choice for the potential U(k), 
i.e. the definition of eq. J2J is extended to all momenta k, 
below and above the Fermi momentum kp. More details 
can be found in review papers or in textbooks [Tlj . We 
extend our calculations up to the three-hole line level of 
approximation, which has been shown to be accurate in 
refs. [2l|, H3] for a variety of NN potentials, including 
the Argonne V14 and vis and the Paris interactions. 

However, since Doleschall's NL potential is used here 
for the first time in BBG calculations of nuclear matter, 
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FIG. 5: Potential energy per particle for symmetric nu- 
clear matter calculated with the non-local potentials NLl 
(squares), which includes the non- locality in the So channel 
only, and NL2 (triangles), which has the non-locality both in 
the 1 So and in the 3 S — D\ channels. Empty symbols cor- 
respond to the two hole-line calculations, while full symbols 
represent the full two + three hole-line results. 



before making comparisons of different EOS, we briefly 
study the accuracy of the BBG expansion in this case. 
This is done in fig. El where the contributions to the 
binding potential energy of the two hole-line (Brueckner) 
and three hole-line diagrams are reported. To keep the 
calculations reliable, the latter has to be substantially 
smaller than the former. 

Two different choices for the NL potential are consid- 
ered: the curves labeled as NLl have been obtained using 
Doleschall's IS potential of ref. only in the 1 Sq chan- 
nel, while for the curves labeled as NL2 the non-local 
potential has been included also in the 3 Si- 3 Di channel. 
In the remaining channels the vig interaction has been 
used. 

The three hole-line contributions turn out to be slightly 
larger than in the case of local potentials (see [2lL 123 . 
[2^). However these corrections always remain within 2- 
3.5% for the NLl case and within 2.5-4% for NL2, thus 
indicating that our calculations up to three hole- lines can 
be considered accurate. It is interesting to note that 
with the non-local potentials the corrections to the BHF 
curves are positive, contrary to the case of the vis, where, 
in the region 1.4 < /cp < 2 fm _1 the three hole-line 
contributions tend to make the EOS more attractive. 



B. Nuclear matter saturation properties 

Having established the validity of our calculations, we 
can now turn our attention to the study of the saturation 
properties of nuclear matter. To better understand the 
role that the different potentials we studied in section ITTI 
play in the nuclear matter EOS it is useful to rewrite the 



Bethe-Goldstone equation @ in the channel c, 

(kl\G c (p,u;)\k'l') = V^(k,k'2 

Q(q) 



(5) 



q 2 dqV^(k,q) ^ (ql"\G c (p,Lo)\k'l') , 
11 LU-E(q) 



where, following standard procedures [Tj|. we have intro- 
duced the angle averaged expressions of the Pauli oper- 
ator, Q, and of the energy, E, in the denominator of the 
integral term. Eq. JSJ closely resembles the Lippmann- 
Schwinger equation the difference being due to the 
presence, in the former, of the Pauli blocking opera- 
tor and of a different energy denominator. Both these 
terms produce a quenching of the integral term contri- 
bution in the Bethe-Goldstone equation with respect to 
the Lippmann-Schwinger case. Therefore the behavior 
of the G-matrix is expected to be rather similar to that 
of the R-matrix, with this quenching effect being weaker 
for the non-local potentials, for which the integral term 
is already small. 

Since the non-locality in the NL potential is introduced 
in the 1 5'o and in the coupled 3 S\- 3 Di channels, we begin 
by studying the contributions of these channels to the po- 
tential binding energy per particle. This is done, at the 
Brueckner level, in fig. El where the results corresponding 
to the non-local CD Bonn (empty dots/dashed line) and 
NL potentials (triangles/dot-dashed) are compared with 
the vis (crosses/solid). In agreement with the discussion 
in section m where we compared bare potentials and R- 
matrix elements, we observe that in the channel the 
CD Bonn and NL potential contributions are very close 
to each other and that they provide additional binding, 
with respect to the vig case. When we consider the 3 Si 
and 3 Di channels, again we see that the non-local po- 
tentials provide more binding, but in this case the NL 
curve is significantly lower than the CD Bonn one. Since 
the contributions from all other channels are very simi- 
lar for the three potentials (with differences of the order 
of 1 MeV at the largest densities considered here), the 
differences observed in the 1 So and 3 Si and 3 D\ chan- 
nels determine the behavior of the total binding potential 
energy, shown in the bottom panel of fig. El 

In order to get further insight on the origin of the dif- 
ferences in the potential binding energy, in fig. [3 we also 
show the potential U(k) for three different choices of the 
density, namely hp — 1.4, 1.6 and 1.8 fm -1 . The most 
important contribution to the total binding energy comes 
from the low k region, where the NL potential is, indeed, 
lower than the others. At large k the behavior of the CD 
Bonn curve differs from the others, which is not surpris- 
ing considering the different high momentum tail of the 
potential V shown, for example, in fig. ^ However this 
behavior does not seem to affect the final result for the 
binding energy. 

Finally, the saturation curve of nuclear matter is shown 
in figurcEl where, in order also to make contact with the 
existing literature, e.g. [T4L ITHl Hfjj. in the upper panel 
we show results at the Brueckner level, while in the lower 
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FIG. 6: Contributions of the different channels to the poten- 
tial binding energy per particle (in MeV) vs. the Fermi mo- 
mentum kF- The upper panel corresponds to the 1 So channel, 
the middle panel to the 3 Si + 3 D\ channels, while the lower 
panel shows the total potential binding energy. All curves are 
calculated at the 2 hole-line level. 



panel we include contributions at the three-hole line level. 

Again, as a reference curve we take the one obtained 
using the Argonne Vis interaction (solid line): as it is 
well established, the saturation density it gives is slightly 
larger than the empirical one, and three-body forces can 
be introduced to correct for this drawback [ll], [24J . This 
is just the issue we want to address, and therefore all the 
calculations will be presented without the introduction of 
three-body forces, in order to see to what extent the non- 
locality can affect this result and in which direction. The 
discrepancy of the calculated saturation point from the 




FIG. 7: Potential U{k) (MeV) of eq. QJ for k F = 1.4 (upper 
panel), 1.6 (middle) and 1.8 fm" 1 (lower) and for the vis 
(solid line), non-local (dot-dashed) and CD Bonn (dashed) 
potentials. 



empirical one can be taken as a measure of the amount 
of three-body forces which is needed. Also relativistic ef- 
fects can shift the saturation point, as described within 
the Dirac-Brueckner (DB) approach Q]. However, it has 
been shown in ref. |25| that the relativistic effects intro- 
duced in DB calculations can be described in terms of 
non-relativistic three-body forces due to virtual creation 
of nucleon-antinucleon pairs. 

Let us now consider the other two interactions, which 
display different degrees of non-locality. The CD Bonn 
at the Brueckner level has been already studied in ref. 
[26j | and more recently in ref. [27|. The corresponding 
saturation curve reported in fig. [S] is in agreement with 
previous studies. The saturation density and energy are 
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FIG. 8: Energy per particle (in MeV) for symmetric nu- 
clear matter calculated for the vis and CD Bonn potentials 
and for the non-local interactions NL1, which includes the 
non-locality in the So channel only, and NL2, where the 
non- locality is included both in the 1 So and in the 3 S — Di 
channels. The upper panel shows results at the two hole line 
level, while in the lower panel three hole-line contributions 
are included. 



substantially larger than for vis. 

The introduction of three hole-line contributions 
slightly improves the situation, since the three hole-line 
diagrams give an overall additional repulsion, which how- 
ever is much smaller than the two hole-line contribution. 
As shown in the bottom panel of fig. 00 the binding en- 
ergy decreases, but still remaining too large with respect 
to the empirical value, while the saturation density is es- 
sentially unchanged. All that indicates that, in order to 
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get the correct saturation point, one needs an amount 
of three-body force larger than in the vi§ case. This 
is somehow at variance with the results for few nucleon 
systems, where the CD Bonn seems to produce binding 
energies closer to the experimental ones [l9|. 

The Doleschall' s potential moves further in this di- 
rection. As in fig[S]we consider again the NL1 and NL2 
cases. The saturation density, as one can see in fig. [SJ is 
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IV. CONCLUSION 

We have analyzed the effects of the non-locality of the 
NN interaction on the saturation curve of nuclear matter 
by considering a set of interactions with different degree 
of non-locality. While the introduction of non-locality is 
able to improve the agreement with phenomenology in 
few-body systems, where it can reduce or even eliminate 
the need of three-body forces, in nuclear matter we ob- 
serve a parallel shift of the saturation point in the wrong 
direction, away from the empirical saturation point. The 
presence of non-locality seems to be unable to solve one of 
the fundamental problem in the many-body description 
of nuclear systems, i.e. the (in)consistency of the three- 
body forces needed in few-body nuclear systems and in 
nuclear matter. 
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